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Quantum signatures of Chimera states 
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Chimera states are complex spatiotemporal patterns in networks of identical oscillators, charac¬ 
terized by the coexistence of synchronized and desynchronized dynamics. Here we propose to extend 
the phenomenon of chimera states to the quantum regime, and uncover intriguing quantum signa¬ 
tures of these states. We calculate the quantum fluctuations about semiclassical trajectories and 
demonstrate that chimera states in the quantum regime can be characterized by bosonic squeezing, 
weighted quantum correlations, and measures of mutual information. Our findings reveal the rela¬ 
tion of chimera states to quantum information theory, and give promising directions for experimental 
realization of chimera states in quantum systems. 
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In classical systems of coupled nonlinear oscillators, 
the phenomenon of chimera states, which describes the 
spontaneous emergence of coexisting synchronized and 
desynchronized dynamics in networks of identical ele¬ 
ments, has recently aroused much interest [1]. These in¬ 
triguing spatio-temporal patterns were originally discov¬ 
ered in models of coupled phase oscillators. In this case, 
they exist due to nonlocal coupling between identical el¬ 
ements of the ensemble [2, 3]. There has been extensive 
work on the theoretical investigation of chimera states [4- 
16], followed by spectacular experimental observations of 
chimera states in optical [17], chemical [18], mechani¬ 
cal [19, 20], electronic [21, 22], optoelectronic [23], and 
electrochemical [24, 25] setups. 

While synchronization of classical oscillators has been 
well studied since the early observations of Huygens in 
the 17th century [26], synchronization in quantum me¬ 
chanics has only very recently become a focus of interest. 
For example, quantum signatures of synchronization in a 
network of globally coupled Van der Pol oscillators have 
been investigated [27, 28]. Related works focus on the 
dynamical phase transitions of a network of nanomechan¬ 
ical oscillators with arbitrary topologies characterized by 
a coordination number [29], and the semiclassical quan¬ 
tization of the Kuramoto model by using path integral 
methods [30]. 

Contrary to classical mechanics, in quantum mechanics 
the notion of phase-space trajectory is not well defined. 
As a consequence, one has to define new measures of 
synchronization for continuous variable systems like op¬ 
tomechanical arrays [29]. These measures are based on 
quadratures of the coupled systems and allow one to ex¬ 
tend the notion of phase synchronization to the quantum 
regime [31]. Additional measures of synchronization open 
intriguing connections to concepts of quantum informa¬ 
tion theory [32], such as decoherence-free subspaces [33], 
quantum discord [34], entanglement [35, 36], and mutual 
information [37]. Despite the intensive theoretical inves¬ 
tigation of quantum signatures of synchronized states, to 
date, studies of the quantum manifestations of chimera 
states are still lacking. 


In this Letter we study the emergence of chimera states 
in a network of coupled quantum Van der Pol oscillators. 
Unlike in previous work [38], we address here the fun¬ 
damental issue of the dynamical properties of chimera 
states in a continuous variable system. Considering the 
chaotic nature of chimera states [11], we study the short- 
time evolution of the quantum fluctuations at the Gaus¬ 
sian level. This approach allows us to use powerful tools 
of quantum information theory to describe the correla¬ 
tions in a nonequilibrium state of the system. We show 
that quantum manifestations of the chimera state ap¬ 
pear in the covariance matrix and are related to bosonic 
squeezing, thus bringing these signatures into the realm 
of observability in trapped ions [27], optomechanical ar¬ 
rays [29], and driven-dissipative Bose-Einstein conden¬ 
sates [39]. We find that the chimera states can be char¬ 
acterized in terms of Renyi quantum mutual informa¬ 
tion. Our results reveal that the mutual information for 
a chimera state lies between the values for synchronized 
and desynchronized states, which extends in a natural 
way the definition of chimera states to quantum mechan¬ 
ics. 

The model .— Similarly to Ref. [27], we consider a 
quantum network consisting of a ring of N coupled Van 
der Pol oscillators. Such a network can be described by 
the master equation for the density matrix p{t) [40] 


p =-Uh , p] + 2'Y^ \K{D{a\) + K2V{a?i) 

1 = 1 


( 1 ) 


where aj , ai are creation and annihilation operators of 
bosonic particles and T>(0) = OpO^ — ^{O^OpP pO^O) 
describes dissipative processes with rates k,i^k ,2 > 0. In 
addition, we have imposed periodic boundary conditions 
ai = for the bosonic operators. In contrast to 

Ref. [27], we consider a nonlocal coupling between the 
oscillators. Therefore, the Hamiltonian in the interac- 
tion picture reads H =% Ya=i 

where V is the coupling strength, d is the coupling range, 
and m 7^ / in the second sum. This kind of coupling im¬ 
plies that Eq. (1) has a rotational symmetry. In the 
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FIG. 1. Space-time representation of the classical chimera 
state for oscillators ai{t) = : (a) 0z(t) and (b) rf{t). 

Parameters: d = 10, n 2 = 0.2avi, V = 1.2a^i, and N = 50. 


particular case d = (A^ + l)/2, N odd, one has all-to-all 
coupling and recovers the results of Ref. [27]. 

Gaussian quantum fluctuations about semiclassical tra¬ 
jectories . — Our aim in this section is to discuss the 
quantum fluctuations about a semiclassical trajectory 
in a way similar to Ref. [41]. We define the expan¬ 
sion of the bosonic operator aflt) = aflt) + d/, where 
a(t) = [(ai(t),..., (aAr(t)] is the semiclassical trajectory, 
and a = (di,..., uat) describes the quantum fluctuations. 
In this work we consider the semiclassical regime, where 
the magnitude of aflt) is much larger than the quantum 
fluctuations described by di. 

Let us consider the displacement operator D [a(t)] = 


exp 


a(t) • a 


• a 


which enables us to define co¬ 


herent states \o:i{t)) = D [cx{t)] \0i), where |0/) is the 
vacuum state of the l-th oscillator, and aflt) \o:i{t)) = 
aflt) \o:i{t)) [42]. By using the expansion of the mas¬ 
ter equation about the mean-field a(t) described in the 
supplemental material [43] , we obtain a master equation 
for the density operator in a co-moving frame p«(t) = 
Dt [a{t)]p{t)D[a{t)] 


N 

^ + 2 ^ [K,V{al) + 4K2\ai\^V{ai) . 

1 = 1 

( 2 ) 

In addition, the coherent dynamics of the fluctuations is 
governed by the Hamiltonian 



1=1 1=1 m=l—d 

m^l 


( 3 ) 

The mean fields appearing in Eq. (3) satisfy the equa¬ 
tion of motion 


l-\-d 

aiit) = ai{t){Ki - 2K2\ai{t)\‘^) - i— am{t) (4) 

m=l—d 

m^l 


with a similar equation for dcl{t). The equations of 
motion Eq. (4) resemble a system of coupled Stuart- 


Landau oscillators [14]. By solving the equations of mo¬ 
tion Eq. (4), one obtains the time-dependent mean field 
exit). Such a mean field plays a fundamental role in the 
description of the master equation Eq. (2). In particular, 
the mean field drives coherent effects such as squeezing 
in Eq. (3) and it determines the time-dependent rates 
which appear in Eq. (2). 

The Classical Chimera state .— Erom our previous 
discussion, the classical equations of motion (4) must be 
satisfied in order to describe the physics in the co-moving 
frame. In the polar representation aflt) = ri{t)e^^^^^^ the 
equations of motion couple amplitude rflt) and phase 
(jflt) of the individual oscillators. We numerically solve 
Eq. (4) for a network of = 50 coupled oscillators 
with coupling range d = 10, considering initial condi¬ 
tions \ai{to)\ ~ ro, where ro = 1.58, and phases drawn 
randomly from a Gaussian distribution in space [43]. Eig- 
ure 1 depicts the time evolution of a classical chimera 
state. In Eig. 1 (a) we show the space-time represen¬ 
tation of the phases (jflt) of the individual oscillators. 
One can observe that for a fixed time, there is a domain 
of synchronized oscillators that coexists with a domain 
of desynchronized motion, which is a typical feature of 
chimera states. Besides the phase, also the amplitude 
exhibits chimera dynamics as we show in Eig. 1 (b). One 
can observe that the width of the synchronized region 
changes with time. Similarly, the center of mass of the 
synchronized region moves randomly along the ring [11]. 
In the case of the uncoupled system with V = 0, the 
individual oscillators exhibit a limit cycle with radius 
ro = which is depicted in the insets of Eig. 2 by 

the green circle. 

Gaussian fluctuations and the Wigner function. — As 
discussed, the classical equations of motion Eq. (4) ex¬ 
hibit a chimera state. By using the knowledge we have 
about the classical trajectory a(t), we can study the 
quantum fluctuations in the co-moving frame by solving 
the master equation Eq. (2). Eor this purpose, we con¬ 
sider the pure coherent state as an initial density matrix 
p{to) = (<^z(^o)|, where \ai{to)\ « 1.58 and 

we choose the phases as in the left panel of Eig. 2. This 
initial condition corresponds to a fixed time to = 3000//i:i 
in Eig. 1. In the co-moving frame, such initial con¬ 
dition reads pocito) = |0z) (0/1* convenience, 

let us write the bosonic operators ai = {qi + ipi)/V^ 
and di = {qi + ip/)/v^ in terms of position and mo¬ 
mentum operators. In terms of complex variables zi = 
(« + ipi)/v^, zi = {qi + ipi)/V^, we define the coor¬ 
dinates = {zi,... ^ zn) in the laboratory frame and 
= {zi,..., zjsf) in the co-moving frame, such that 
z = a(t) + z. The variables qi^qi and pi^pi denote posi¬ 
tion and conjugate momentum, respectively. 

We consider the Wigner representation Woc{Fi^ t) of the 
density operator p(x{t)^ where = (gi,pi,..., gA^ 5 PA^)• 
By using standard techniques of quantum optics [40], 
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FIG. 2. Quantum signatures of the classical chimera state, 
(a) Snapshot of the phase chimera depicted in Fig. 1 at 
= 3000. We consider an initial density matrix p(to) 
which is a tensor product of coherent states centered around 
the positions of the individual oscillators as depicted in the 
insets (Husimi function), (b) After a short-time interval 
KiAt = 0.5, quantum correlations appear in the form of 
squeezing (black double arrows in the insets). Parameters: 
d = 10, = 0.2/^i, V = 1.2/^i, and N = 50. 


the master equation Eq. (2) can be represented as a 
Fokker-Planck equation for the Wigner function which 
depends on the mean field solution of Eq. (4) and con¬ 
tains information of the chimera state. In the supple¬ 
mental material [43], we provide the explicit form of 
the Eokker-Planck equation for IFq-. Eortunately, even 
though the coefficients of the equation are time de¬ 
pendent, one can derive an exact solution Wcc{R,t) = 

(27r)“^(detexp , where ^{t) 

is the covariance matrix, whose matrix elements ^ij = 
(^^{RiRj ^ RjRi)'^ — informa¬ 

tion about the correlations between quantum fluctua¬ 
tions R 21-1 = Qi and R 21 = pi. The angular brackets 
= tr(pQ,0) denote the expectation value of an op¬ 
erator O calculated with the density matrix po,. 

The solution Woc{RR) corresponds to a Gaussian dis¬ 
tribution centered at the origin in the co-moving frame. 
In the laboratory frame, the Wigner function is centered 
at the classical trajectory cx{t). However, due to the 
chaotic nature of the classical chimera state [11], our ex¬ 
act solution is just valid for short-time evolution. The 
Husimi function Q{z) = ^{z\p{t)\z) is intimately re¬ 
lated to the Wigner function [40, 43] and can be ob¬ 
tained numerically by using the Gutzwiller ansatz [27]. 
The insets in the right panel of Eig. 2 depict the Husimi 
functions of the individual nodes after a short evolution 
time At = 0.5//^i. One can observe that even if one 
prepares the system in a separable state, quantum fluc¬ 
tuations arise in the form of bosonic squeezing of the os¬ 
cillators [40]. In the insets of Eig. 2, the arrows indicate 


the direction perpendicular to the squeezing direction for 
the individual oscillators. Eor oscillators within the syn¬ 
chronized region, the squeezing occurs almost in the same 
direction. In contrast, the direction of squeezing is ran¬ 
dom for oscillators in the desynchronized region, which 
reflects the nature of the chimera state. 

Quantum signatures of a chimera state in the covari¬ 
ance matrix .— Now let us study the consequences of 
the exact solution for the short-time evolution of the 
Wigner function. Once we obtain the solution of the 
equations of motion Eq. (4) , one can find the correspond¬ 
ing covariance matrix ^{t). As we have defined in the 
introduction, a chimera state is characterized by the co¬ 
existence in space of synchronized and desynchronized 
motion. Therefore, to understand the quantum mani¬ 
festations of a chimera state, we require to study also 
quantum signatures of synchronized and desynchronized 
dynamics. 

The left column of Eig. 3 show snapshots of the phases 
for (a) chimera, (b) synchronized, and (c) desynchro¬ 
nized mean-field solutions of Eq. (4). The central col¬ 
umn of Eig. 3 depicts the corresponding covariance ma¬ 
trices after a short evolution time At = 0.5//^i. Eor 
every plot, we have initialized the system at time ti as 
a tensor product of coherent states \ai{ti)) centered at 
the positions aptQ of the individual oscillators. As a 
consequence, the covariance matrix at the initial time is 

diagonal ^ 2 i-i, 2 i-i{ti) = = ^/2 and ^ 2 Z, 2 z(ti) = 

(p‘i'^ = ^/2, which reflects the Heisenberg uncertainty 

principle because = (pi"^ = 0. 

After a short evolution time, quantum correlations are 
built up due to the coupling between the oscillators, 
and the covariance matrix exhibits a nontrivial struc¬ 
ture which is influenced by the mean field solution. Eor 
example, the central panel of Eig. 3 a) shows a ma¬ 
trix plot of the covariance matrix corresponding to a 
chimera state obtained from the same initial condition 
as in Eig. (2). The covariance matrix acquires a block 
structure, where the upper 40 x 40 block (corresponding 
to nodes / = !,..., 20) shows a regular pattern matching 
the synchronized region of the chimera state. Similarly, 
the lower 60 x 60 block shows an irregular structure which 
corresponds to the desynchronized dynamics of the oscil¬ 
lators / = 21,...,50. In a similar fashion, Eigs. 3 b) 
and c) show the matrix ^ for completely synchronized 
and desynchronized states, respectively. In the case of 
a chimera state, this coincides with the results shown in 
Eig. 2, where the squeezing direction of the oscillators 
is related to the classical solution. In order to quantify 
these observations we define the weighted correlation as 

XT l-\-d 

= ^ E (5) 

m=l—d 

m^l 
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FIG. 3. Quantum fluctuations after a short-time evolution. 
Similarly to Fig. 2, we consider an initial density matrix 
p{ti) which is a tensor product of coherent states centered 
around the classical positions of the oscillators. Snapshots of 
the phase (left column) and covariance matrices (central col¬ 
umn) after short-time evolution At = 0.5 of the states: (a) 
chimera for V = (b) synchronized state for V = 

and (c) desynchronized state for V = 0.8/^i. Right column: 
Weighted spatial average ^i(t) of the covariance matrix for 
the states shown in a), b) and c), respectively. Parameters 
d = 10, k ,2 = 0.2/^i, and N = 50. 


This spatial average highlights the structure of the co- 
variance matrix. The right column of Fig. 3 shows T/(t) 
for a) chimera, b) synchronized, and c) desynchronized 
states. The chimera state exhibits a regular and an ir¬ 
regular domain, exactly as the classical chimera does. 

Chimera states and Renyi quantum mutual 
information .— Now let us consider a partition of 
the network into spatial domains of size L and N — 
which we call Aliee (A) and Bob (B), respectively. 
This partition can be represented by considering a 
decomposition of the covariance matrix 




^A{t) ^AB{t) A 
^AB (^) {t) ) 


(6) 


To study the interplay between synchronized and desyn¬ 
chronized dynamics, which is characteristic of a chimera 
state, we propose the use of an entropy measure [32, 37]. 
Of particular interest is the Renyi entropy S^{p) = 
(1 — /i)“^ Intr(p^), /i G N, of the density matrix p, which 
is discussed in Ref. [44]. In terms of the Wigner rep¬ 
resentation of poc^ the Renyi entropy for /i = 2 reads 


S 2 {poc) = - In 


fW^(R,t)d^^R 


Now let us consider 


the bipartite Gaussian state pab = Pa composed of Aliee 
and Boh and define the tensor product ppef = Pa O Pb of 
the two marginals pA and pe- 

To measure Gaussian Renyi-2 mutual information 
2^2(M:b) = S2{pa) + S2{pb) “ 5 ' 2 (/ 9 ab), we require the 
calculation of the relative sampling entropy between the 



FIG. 4. Renyi quantum mutual information for the states 
shown in Fig. 3. The green dots, blue diamonds, and purple 
triangles represent the chimera, synchronized, and desynchro¬ 
nized states, respectively, a) Gaussian Renyi-2 mutual infor¬ 
mation X 2 (pa:b) as a function of the size L of Alice after an 
evolution time At = 0.5//^i. b) The time evolution of the mu¬ 
tual information during the time interval At for a fixed size 
Lc — 20. Inset: scheme of the nonlocally coupled network. 
Parameters: d = 10, k ,2 = 0.2/^i, and N = 50. 


total density matrix pab and the reference state ppef as 
shown in Ref. [44]. This leads to a formula X 2 (pa:b) = 
I In (det det / det in terms of the covariance ma¬ 
trix Eq. (6). Figure 4 a) shows the variation of X 2 (pa:b) 
as a function of the size L of the partition after an evo¬ 
lution time At = 0.5/ki. One can observe that for a 
chimera state the mutual information is asymmetric as a 
function of L and there is a critical size Lc = 20, where 
a dramatic change of the correlations occurs. 

Now let us consider the chimera state shown in Fig. 2, 
and let us consider a partition where the size of Aliee is 
Lc = 20. Fig. 4 b) shows the time evolution of mutual 
information for such a state. In addition, by using the 
same partition as for the chimera state, we calculate the 
mutual information for the synchronized and desynchro¬ 
nized states depicted in Figs. 3 b) and c), respectively. 
Our results reveal that the chimera state has a mutual in¬ 
formation which lies between the values for synchronized 
and desynchronized states. This resembles the definition 
of a chimera state given at the beginning of the article. 

Conelusion .— We have shown that quantum signa¬ 
tures of chimera states appear in the squeezing of coher¬ 
ent states, in the covariance matrix, and in measures of 
mutual information. To quantify the structure of the co- 
variance matrix, we have introduced a spatial average of 
the quantum correlation, which reveals the nature of the 
classical trajectory, i.e., chimera, synchronized, or desyn¬ 
chronized state. The mutual information for a bipar¬ 
tite state X 2 (pa:b) extends the definition of a chimera to 
the quantum regime and highlights the relation to quan¬ 
tum information theory. A possible experimental real¬ 
ization of our model could be carried out by means of 
trapped ions [45], as it was suggested in Ref. [27]. Other 
experimental possibilities include Bose-Einstein conden¬ 
sation in the presence of dissipation and external driv¬ 
ing [39, 46]. In this context our approach is particularly 
interesting, because the continuum limit of the mean field 
Eq. (4) is a complex Ginzburg-Landau equation, which 
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is nonlocal in space [2]. In this sense, our linearized mas¬ 
ter equation Eq. (2) enables us to study the Bogoliubov 
excitations around the mean field solution. 

Acknowledgments .— V.M. Bastidas thanks L. M. Va¬ 
lencia and Y. Sato. The authors acknowledge inspiring 
discussions with J. Cerrillo, S. Restrepo, G. Schaller, and 
P. Strasberg. This work was supported by DFG in the 
framework of SFB 910. 


* victor@physik.tu-berlin.de 

[1] M. J. Panaggio and D. M. Abrams, Nonlinearity 28, R67 
(2015). 

[2] Y. Kuramoto and D. Battogtokh, Nonlin. Phenom. in 
Complex Syst. 5, 380 (2002). 

[3] D. M. Abrams and S. H. Strogatz, Phys. Rev. Lett. 93, 
174102 (2004). 

[4] G. C. Sethia, A. Sen, and F. M. Atay, Phys. Rev. Lett. 
100, 144102 (2008). 

[5] C. R. Laing, Phys. D 238, 1569 (2009). 

[6] A. E. Hotter, Nat. Phys. 6, 164 (2010). 

[7] E. A. Martens, C. R. Laing, and S. H. Strogatz, Phys. 
Rev. Lett. 104, 044101 (2010). 

[81 S. Olmi, A. Politi, and A. Torcini, Europhys. Lett. 92, 
60007 (2010). 

[9] G. Bordyugov, A. S. Pikovsky, and M. G. Rosenblum, 
Phys. Rev. E 82, 035205 (2010). 

[10] J. H. Sheeba, V. K. Chandrasekar, and M. Lakshmanan, 
Phys. Rev. E 81, 046203 (2010). 

[11] O. E. Omel’chenko, M. Wolfram, and Y. Maistrenko, 
Phys. Rev. E 81, 065201(R) (2010); M. Wolfram and O. 
E. Omel’chenko, Phys. Rev. E 84, 015201 (R) (2011). 

[12] 1. Omelchenko, Y. L. Maistrenko, P. Hovel, and E. Scholl, 
Phys. Rev. Lett. 106, 234102 (2011). 

[13] 1. Omelchenko, O. E. Omel’chenko, P. Hovel, and E. 
Scholl, Phys. Rev. Lett. 110 , 224101 (2013). 

[14] A. Zakharova, M. Kapeller and E. Scholl, Phys. Rev. 
Lett. 112 , 154101 (2014). 

[15] P. Ashwin and O. Burylko, Chaos 25, 013106 (2015). 

[16] 1. Omelchenko, A. Zakharova, P. Hovel, J. Siebert, and 
E. Scholl, arXiv: 1503.03377 (2015). 

[17] A. M. Hagerstrom, T. E. Murphy, R. Roy, P. Hovel, 1. 
Omelchenko, and E. Scholl, Nat. Phys. 8, 658 (2012). 

[18] M. R. Tinsley, S. Nkomo, and K. Showalter, Nat. Phys. 
8, 662 (2012). 

[19] E. A. Martens, S. Thutupalli, A. Fourriere, and O. Hal- 
latschek, Proc. Natl. Acad. Sci. U.S.A. 110 , 10563 (2013). 

[20] T. Kapitaniak, P. Kuzma, J. Wojewoda, K. Czolczynski, 
and Y. L. Maistrenko, Scientific Reports 4, 6379 (2014). 

[21] L. Larger, B. Penkovsky, and Y. L. Maistrenko, Phys. 
Rev. Lett. Ill, 054103 (2013). 

[22] D. P. Rosin, D. Rontani, N. D. Haynes, E. Scholl, and D. 
J. Gauthier, Phys. Rev. E 90 90, 030902(R) (2014). 

[23] L. Larger, B. Penkovsky, and Y. Maistrenko, 
arXiv: 1411.4483 (2014). 

[24] M. Wickramasinghe and 1. Z. Kiss, PLoS ONE 8, e80586 
(2013). 

[25] L. Schmidt, K. Schonleber, K. Krischer, and V. Garcia- 
Morales, Chaos 24, 013102 (2014). 

[26] M. Kapitaniak, K. Czolczynski, P. Perlikowski, A. Ste- 
fanski, and T. Kapitaniak, Phys. Rep. 517, 1 (2012). 


[27] T. E. Lee, and H. R. Sadeghpour, Phys. Rev. Lett. Ill, 
234101 (2013). 

[28] S. Walter, A. Nunnenkamp, and C. Bruder, Phys. Rev. 
Lett. 112 , 094102 (2014). 

[29] M. Ludwig and F. Marquardt, Phys. Rev. Lett. Ill, 
073603 (2013). 

[30] 1. H. de Mendoza, L. A. Pachon, J. Gomez-Gardehes, and 
D. Zueco, Phys. Rev. E 90, 052904 (2014). 

[31] A. Mari, A. Farace, N. Didier, V. Giovannetti, and R. 
Fazio, Phys. Rev. Lett. Ill, 103605 (2013). 

[32] J. Eisert and M.B. Plenio, Int. J. Quantum. Inform. 
01, 479 (2003); S. L. Braunstein and P. van Loock, Rev. 
Mod. Phys. 77, 513 (2005); C. Weedbrook, S. Pirandola, 

R. Garcfa-Patron, N. Cerf, T. Ralph, J. Shapiro, and S. 
Lloyd, Rev. Mod. Phys. 84, 621 (2012). 

[33] G. Manzano, F. Galve, G. L. Giorgi, E. Hernandez- 
Garcfa, and R. Zambrini, Sci. Rep. 3, 1439 (2013). 

[34] G. Manzano, F. Galve, and R. Zambrini, Phys. Rev. A 
87, 032114 (2013). 

[35] T. E. Lee, C.-K. Chan, and S. Wang Phys. Rev. E 89, 
022913 (2014). 

[36] G. L. Giorgi, F. Plastina, G. Francica, and R. Zambrini, 
Phys. Rev. A 88, 042115 (2013). 

[37] V. Ameri, M. Eghbali-Arani, A. Mari, A. Farace, F. 
Kheirandish, V. Giovannetti, and R. Fazio, Phys. Rev. 
A 91, 012301 (2015). 

[38] D. Viennot, L. Aubourg, arXiv: 1408.4585 (2014). 

[39] L. M. Sieberer, S. D. Huber, E. Altman, and S. Diehl, 
Phys. Rev. Lett. 110 , 195301 (2013); U. C. Tauber and 

S. Diehl, Phys. Rev. X 4, 021010 (2014) . 

[40] H. J. Carmichael,, Statistical Methods in Quantum Optics 
2 (Springer, Berlin, 2007). 

[41] D. A. Rodrigues and A. D. Armour, Phys. Rev. Lett. 104, 
053601 (2010); N. Lorch, J. Qian, A. Clerk, F. Marquardt, 
and K. Hammerer, Phys. Rev. X 4, 011015 (2014). 

[42] K. E. Cahill and R. J. Glauber, Phys. Rev. 177, 1857 
(1969). 

[43] See supplemental material for additional details. 

[44] G. Adesso, D. Girolami, and A. Serafini, Phys. Rev. Lett. 
109, 190502 (2012). 

[45] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. 
Mod. Phys. 75, 281 (2003). 

[46] J. Kasprzak, M. Richard, S. Kundermann, A. Baas, P. 
Jeambrun, J. M. J. Keeling, F. M. Marchetti, M. H. 
Szymahska, R. Andre, L. Staehli, V. Savona, P. B. Lit- 
tlewood, B. Deveaud, and Le Si Dang, Nature (London) 
443, 409 (2006). 



6 


SUPPLEMENTAL MATERIAL FOR “QUANTUM SIGNATURES OF CHIMERA STATES” 

V. M. Bastidas^, 1. Omelchenko^, A. Zakharova^, E. Scholl^, and T. Brandes^ 

^Institut fiir Theoretische Physik, Technische Universitat Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 


CALCULATION OF THE GAUSSIAN QUANTUM FLUCTUATIONS FOR A SINGLE VAN DER POL 

OSCILLATOR 


In the case of a single Van der Pol oscillator [1, 2] we can write the master equation Eq. (1) as follows 

p{t) = -^[Hint,p] + 2 ki (^a^pa- ^{p,aa^}^ +2 k 2 (^a'^p{a^f - , ( 7 ) 

where in the main text we set Hint = 0 for the single oscillator and only consider coherent dynamics through the 
coupling. 

Eollowing a similar method as the one used in Ref. [3], if one defines the reduced density matrix in the co-moving 
frame pa{t) = & [<^(t)] p{t)D we obtain a master equation with Liouville operators £i and £2 


Pa{t) = -^[H^°‘\t),pa{t)] +D^ [a{t)] 2ki (^d^pa- ^{p,dd^}^ + 2k2 (^a^p{d^f - ]^{p,{dtfa^} 

= — pcx{t)] ^ Clpcx ^ C2PCX 1 


D[a{t)] 


( 8 ) 


where we have defined [<^(t)] (^int ~ ihdt)D [o:{t)] and the anticommutator {A^ B} = AB + BA. Eor 

later purposes, we need to use the identity 


[a{t)] dtD [ct{t)] = ^[d(t)a*(t) — (a(t)d*(t)] + ih[a{t)a^ — d*(t)a] 


(9) 


We are interested in the action of the displacement operator D [a{t)] on the other terms of the master equation. 
By using elementary properties of the displacement operator we obtain 


jCipa{t) = [a{t)] (d^pd - ^{p, aa'*'} ) D [a{t)] 


= 2/^1 (dVad- ^ [ihKiad\pa] - ^ [-ihKia^d, pa] • 


( 10 ) 


However, the dissipative term which is proportional to /^2 requires more attention. Therefore, we study each term 
individually 


C2Pa{t) = 2 k2D^ [a(i)] (d^p{dtf - b/9, (at)2aq) 1) [a(t)] 


2k2 + 4k2Q;* - ]j^{pa, + 4K2a (^pa{dtf‘ - 


+ 8k2|q;P [dpad'^ - hpaja^a} ) - f [ihK2{a*)‘^d^, pa] - f [-\fiK,2a^{d'^)‘^ , Pa] 


— — [2ihK2a{a*)'^d, Pa] — — [—2ihK2a*a‘^d^,pa] 


( 11 ) 


Now we compare with the general form of the Lindblad master equation [4] with Lindblad operators 

Kt) = -^[-f7int,p] + (l^pLI - ^{p,lIlA , 

II ^ ' 


( 12 ) 
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which enables us to introduce an effective Hamiltonian to describe the dynamics between quantum jumps 

fl 

In order to study quantum fluctuations about semiclassical trajectories, we consider an effective non-hermitian Hamil¬ 
tonian as an expansion in powers of quantum fluctuations 

2i/^2<^*d^d^ — 2i/^2<^(d^)^d) 
h{—iKiaa^ — 4i/i:2|<^pd^d + i/T:2(<^*)^d^ — i/T:2<^^(d^)^) 

h{—i[a{t)a^ — a*{t)a] + ii^iaa^ — i^ia^'a + 2i/i:2Q^(<^*)^d — 2i/i:2<^*Q^^d^) • (14) 

In terms of the effective Hamiltonian, we can write the master equation Eq. (8) as follows 

Pa{t) = -^H^'gpait) + ^pa{t){Hl^)* + 2 k 20? Paio^f' +‘iK2a* p^o) + ‘iK2aapa{a^)‘^ 

+ 2/^1 Peth + 8/^21I . (15) 


Now we have all the necessary elements in order to study the quantum fluctuations. First, we eliminate the linear 
terms in the Hamiltonian. This is achieved as long as the condition 

a{t) = — 2i<i2a{t)\a{t)\^ (16) 


is satisfied. Eq.(16) corresponds to the equation of motion of a Van der Pol oscillator. In order to perform a 
semiclassical treatment of the master equation [3, 4], we observe that \o:{t)\ ^ 1. This is equivalent to saying that 
the action of the semiclassical system is greater than the action quanta, which allows us to perform a semiclassical 
approximation to the master equation up to quadratic order in the quantum fluctuations 


Pa{t) « pa{t)] +2ki ^a^paCL- ^{pa^ao^}^ +8K2|a|^ {apaO^ - ]^{pa,o)a}^ , (17) 


where 

H^\t) = \hl<i2[oL^ — \hl<i2[oL{t)]^{d)Y • 


(18) 


INITIAL CONDITIONS 

In the case of the uncoupled system V = 0, a single Van der Pol oscillator [1, 2] exhibits a limit cycle with radius 
ro = 'sj = 1-58 for the parameters /^2 = 0.2/^i used in the main paper. For convenience, in the coupled system 
we consider initial conditions at t = 0 in such a way that each oscillator has the same amplitude |(az(0)| ~ 1.58. In 
addition, we consider phases obtained from a Gaussian distribution 0/(0) = exp[—where —247r < 0 < 247r 
is a random number, fi = N/2 and a = 9. Fig. 5a) shows the initial conditions. In terms of the coordinates 
ai{t) = initial conditions must satisfy + P^/^(0) ~ 2.24, which defines the green circle in 

Fig. 5b). 

In the main text, although we consider different coupling strengths V, we use the same initial conditions as in 
Fig. 5 to obtain the chimera, synchronized and desynchronized states. In order to obtain the snapshot of the chimera 
state depicted in Fig. 3a), we let the system evolve up to a time /^ito = 3000.5 for a coupling strength V = 1.2. 
Correspondingly, to obtain the snapshot of the synchronized solution shown in Fig. 3b), we let the system evolve 
a time /^itsyn = 25.5 for V = 1.6. Finally, the snapshot of the desynchronized in Fig. 3c) is obtained after a time 
evolution /^itdesyn = 8000.5 for V = 0.8. 


EXPLICIT FORM OF THE FOKKER-PLANCK EQUATION 


Now let us define the Wigner representation of the density operator p«(t) [4], as follows 






r2N 


g-A-2*+A*-2^^ 


pa{t)e 


— A'Ct ~t“A 'd 


(19) 









TT 



FIG. 5. Initial conditions used in the main paper, a) Initial distribution of the phases 0z(O) drawn randomly from a Gaussian 
distribution in space, b) Phase-space representation of the initial conditions for the oscillators. The green circle represents the 
limit cycle with radius |<az(0)| 1.58, where ai{t) = . Parameters 1, d = 10, k ,2 = 0.2/^i, and iV = 50. 


where A = (Ai,...,A 7 v) denote the integration variables, z = {zi^..., z^) and zi = {qi In addition, 

the variables qi and pi denote position and momentum, respectively. The Wigner function is related to the Husimi 
function in the co-moving frame Qcx{z) = ^ {z \ pcx{t) \z) via the transformation [4] 

Qociz) = ^ j . ( 20 ) 

By using the Wigner representation, the master equation Eq. (2) can be mapped into a Fokker-Planck equation [4] 


dW^ 

dt 


E + {4:K2\ap - Ki)di^zi + (2K2\ap + y) 5? 


1 = 1 


Wo. 


V 


N l+d 




( 21 ) 


1=1 m=l—d 


By defining = (gi,Pi,. •., the Fokker-Planck equation can be written also in terms of the quadratures qi 

and Pi 


f)W 1 


dt 

i=l i=l 

Although this Fokker-Planck equation has time dependent coefficients, one can derive an exact solution [4] 


( 22 ) 


Wc.{R,t) 




(27r)^\/det ^ 

where the covariance matrix ^{t) is a solution of the differential equation ^{t) = + ^(t). 


(23) 
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